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Abstract 

In this paper we study the description of saturation in Balitsky, 
Jalilian-Marian, Iancu, McLerran, Weigert, Leonidov and Kovner (B- 
JIMWLK) picture when restricted to observables made up only from 
dipole operators. We derive a functional form of the evolution equa- 
tion for the dipole probability distribution and find a one-parameter 
family of exact solutions to the dipole evolution equations. 



1 Introduction 

One of the outstanding issues in perturbative QCD is the understanding of 
parton saturation. Standard evolution equations like BFKL lead to ris- 
ing cross sections (or structure functions) which violate the unitarity bound. 
These are modified, when one takes into account recombination effects, into 
an infinite set of coupled nonlinear equations for correlation functions of 
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Wilson line operators j2]. This set of equations is expected to be equiva- 
lent to the Jalilian-Marian, Iancu, McLerran, Weigert, Leonidov and Kovner 
(JIMWLK) functional equations [3], which can be recast in the form of 
Fokker-Planck equations for the (functional) probability distribution of Wil- 
son lines Zy[{f/ X }]. The probability distribution Zy[{£7 x }] encodes expecta- 
tion values of observables made up of the longitudinal Wilson lines U x , where 
a; is a coordinate in the transverse space: 

(f[{U x }\) = J VU X f[{U x }} Z Y [{U X }} (1) 

The full knowledge of Zy[{U x }] is of course nonperturbative, but the evolu- 
tion in rapidity is governed by the JIMWLK functional equation: 

^Z Y [{U X }\ = a.y<?x<?v (iV a Ux xf y iV h Uy ) Z Y [{U X }] (2) 

where the derivatives are covariant functional derivatives w.r.t. the 
unitary matrices while the kernel Xxy ls given by the following expression: 

«t 4/^ (!^f i 1 - - wr) (3) 

See j3] and [5| for further details. Despite the apparent compactness of 
the above equation which summarizes the whole nonlinear infinite set of 
evolution equations for correlation functions it has proved itself to be difficult 
to analyze. Apart from some general properties not much is explicitly known 
about the details of saturation/unitarization in this picture. 

It is interesting, therefore, to study simplified versions. One noteworthy 
example is the Balitsky-Kovchegov (BK) equation [21 El obtained when cor- 
relation functions of traces of products Wilson lines factorize into products of 
expectation values. This leads to a nonlinear equation for the dipole density 

w {Sm) = &I d2%2 £k [{So2) {S2i) " Nc {Soi)] (4) 

where the «Sy are (the S-matrix versions of) the dipole operators 

<% • (5) 

The Balitsky-Kovchegov equation was studied both numerically and analyti- 
cally [7j and boasts universal properties of geometrical scaling [S] (appearing 
here as universal traveling wave solutions jU]). 
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In order to better understand the full hierarchy of JIMWLK equations, it 
is of particular interest to find some intermediate approximations which are 
not so drastic as BK, yet remain simpler than the full JIMWLK hierarchy. 
One approach might be to restrict oneself just to observables made up from 
dipole operators, as these have a very clear intuitive physical meaning, yet to 
refrain from making any other simplifying assumptions. The generalization 
with respect to BK means that one has to allow for nontrivial correlations 
between dipoles i.e. allowing for (S01S23) 7^ (<Soi) (£23) etc. In [TU] we stud- 
ied the first step beyond BK with nontrivial 2-dipole correlations. Other 
approaches going beyond BK include [TTJ [T2J EH . 

The aim of this paper is to derive an analogue of (J2J) restricted to prob- 
ability distributions made up of dipole operators only, and to find a one- 
parameter class of exact solutions of these dipole evolution equations with 
some specific multi-dipole correlations. 

The comparison of the dipole evolution hamiltonian to the full JIMWLK 
one may be a first step in understanding the role of higher multipole operators 
in perturbative unitarization/saturation processes. 

The plan of this paper is as follows. In section 2 we derive the evolution 
hamiltonian for the JIMWLK hierarchy restricted to dipole operators, then 
in section 3 we discuss some of its properties. In section 4 we construct a 
one-parameter solution of the dipole evolution equations and we close the 
paper with a summary and discussion. 



2 The JIMWLK hierarchy restricted to dipoles 

In order to calculate the evolution of an observable of the form {S01S23 . . . <S n _i n ) 
it is most convenient to use Balitsky's approach where the evolution with ra- 
pidity is obtained using contractions between pairs of the unitary matrices 
U x s appearing in the dipole operators and virtual corrections for each U x ac- 
cording to eqs. (119), (120) in |2|. One can verify that contractions between 
U x 's belonging to different dipoles always give rise to multipoint traces, so if 
one wants to keep only terms involving dipoles one performs contractions only 
within each dipole and therefore each dipole operator evolves independently 
through the operator form of the BK kernel: 

d g 2 r (5 \ 

gyfiSxy] = ^3 / d2xd2 y d2z \ Js~f[ Sx y\) K xzy(S xz S zy - N c S xy ) (6) 
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Introducing the dipole probability distribution Zy[{iS xy }] which encodes ex- 
pectation values through 

(f[{S xy }}) = J VS xy f[{S xy }] Z Y [{S xy }] (7) 

we can rewrite the previous equation as 

'K XZ y(S XZ S Z y — N c S X y\Zy [{_S X y}\ (8) 

We will now recast the right hand side as 

J VS xy f[S xy }HZ Y [{S xy }} (9) 

Then the evolution equation for the dipole probability distribution will have 
the form 

^Zy[{S xy }]=HZ Y [{S xy }] (10) 

If we perform a functional integration by parts in (JBJ) and assume the van- 
ishing of boundary terms we obtain an explicit and quite simple form of the 
evolution hamiltonian 7i: 

-gyZy[{S xy }] = — - J d 2 xd 2 yd 2 zj—- {K xzy (S xz S zy - N c S xy )} Z Y [{S xy }] 

xy 

(11) 



3 Some remarks on the dipole hamiltonian 



As an example of how the functional evolution equation (jTTj) works in practice 
let us quickly rederive the equation for 2-dipole correlations obtained in ^0] • 
Let us start from the identity 



o f o 

7^7 (<Soi<S>23) = J VS xy S 01 S 2 3 -g^Z Y [{S xy \] (12) 
Now we use (jlip and perform integration by parts to get 



J VS xy J d 2 xd 2 yd 2 zj^— (S01S23) {K xzy {S xz S zy - N c S xy )} Z Y [{S xy }} 



in 3 



(13) 
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After performing the functional differentiation and using the definition of 
Zy[{i5>a;j,}] (J7|) we get the equation obtained in [TUj : 

d f 

— (iS 02 5 2 'i) = J d 2 x 3 [(S 03 S 32 S2'i) - N c (SmSxi)] K 032 + 

+ [(^^^kSqs) - N c (5 02 5 2 a)] K lwf . (14) 

Let us now contrast some properties of the dipole evolution equation (jTTjl 
with the full JIMWLK one (J2J). The latter one is second order in functional 
derivatives and, as shown by Weigert, the hamiltonian is positive definite, 
with a unique attractive fixed point (Zy[{f/ Z }] = 1). On the other hand 
the dipole equation is first order and as such does not seem to have such 
positivity and uniqueness properties. Even more so, from the form of (|TT|) 
one would expect that an infinite set of zero- modes would exist. This would 
lead to many probability distributions determined by 

nz Y [{s xy }} = o (is) 

which are invariant under evolution. 

We cannot a-priori rule this out but let us show on a simple toy model 
how one can nevertheless escape this conclusion. Assume the following simple 
evolution equation for a dipole operator 

W s = s-s 2 (16) 

which exhibits saturation behaviour quite similar to BK. This leads to the 
hamiltonian for the probability distribution Zy[s]: 



^-Z Y [s\ =d s 
dY 



-s)Z Y [s}\ (17) 
One can now find explicitly the dangerous stationary solutions: 



However these solutions are nonphysical as they are nonnormalizable due to 
the singularities at s — 0, 1, and therefore cannot represent any probability 
distribution. It would be very interesting (but rather difficult) to carry out 
a similar analysis for the real dipole equation (jll|) . 
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4 A family of factorized solutions to the dipole 
hierarchy of evolution equations 



It is interesting to look for exact solutions of the full set of dipole evolution 
equations. It turns out to be easier to find a solution of the type considered 
here using directly a set of evolution equations for the moments of dipole 
operators rather than finding first the probability distribution Zy[{«S a;2/ }]. 

The JIMWLK-Balitsky hierarchy of evolution equations, when restricted 
to dipole operators has the following form: 

d n 

-g^<Si<S2 ■ ■ S n = Si . . . Si-i(7iSi)Si + i . . . S n (19) 

In terms of dipole density operators Hi related to the Si through 

S i = N c {\-M i ) (20) 

the hamiltonian acts as 

HMi = HN lV = ±— ^ / d 2 xK ixV [N ix + N xi > - M w - M ix M xi ] (21) 

Let us look for solutions to this set of coupled equations in a factorized 
form: 

(MM . . . M n ) = c n (M) (M> • • • (AQ = c n N 1 N 2 ...N n (22) 

with c n some coefficients (ci = 1). Note that (. . .) is the ordinary expectation 
value and not just the connected part. If all c n would be equal to one, 
this would imply no nontrivial correlations and would be equivalent to the 
standard BK equation. 

Let us take the expectation value in the evolution equation (fT9"j) : 

— ((l - M) • • • (i - AQ> = - - M) ••• (i - - M + i) 

...(l-K)) (23) 

We may now use first the assumption of factorization (}2*2*|) in the left hand 
side of (|23p. Thus the F-derivative acts on products of single dipole densities 
Ni = (Afi). Then for each iVj we use the evolution operator 

-qY N i = yj J d xK i*i> i N ™ + N *v ~ N a' ~ C2N ix N xt ,] (24) 
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Note the appearance of c 2 w.r.t. (|2T|l as we are taking expectation values. We 
can now insert the operator form of (|21j) into the right hand side of (|23|) and 
then perform the averaging using the factorization ansatz ([22)1 . Comparing 
terms we see that all terms match if the following holds: 

Cj+i = c 2 Ci (25) 

which has a simple solution 1 

c n = cr 1 (26) 

Therefore all multidipole correlation functions reduce to one-dipole densities 
but with an enhancement factor: 

(MM ...AO = q^NtNt ...N n (27) 
And all the one dipole densities N satisfy a modified BK equation 
d a 2 N r 

g^ N n' = J d2 ^i' Wix + N xi , - N u > - c 2 N ix N xl ,} (28) 

with c 2 being an a-priori arbitrary parameter. However in order to have a 
solution which saturates at Ni < 1, c 2 has to be greater or equal to 1. Indeed 
let us note that (|2"7|) with N = l/c 2 is a fixed point of the dipole evolution 
equations. It would be interesting to recast this solution in the language of 
probability distributions Zy \{S xy \\ and analyze its normalizability proper- 
ties. 

We note that the same modified equation was proposed in [TT] to acccount 
for effects of correlations in nuclei on saturation. It is interesting to note 
that the same equation arises from an exact solution of the dipole evolution 
hierarchy. 

Let us note how does the factorized solution of ^U] fit into the above 
general framework. That solution was obtained assuming that the dipole 
evolution equations close at the level of two equations for (Af\) and (MM) 
i.e. that there are no nontrivial 3-dipole correlations. In order to recover 
that solution in our framework we have to analyze the connected correlation 
functions 

(MM • • .M n ) c = d n NN 2 ...N n (29) 

1 After these results were obtained M. Lublinsky informed me that the same factorized 
solution appeared in a framework of generalized BK in ^T] . 
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The coefficients d n are related to the c n 's through 



c 2 = c\ + d 2 

c 3 = c\ + 3cid 2 + d 3 

C4 = c\ + 6c? 2 c^ + 3d 2 + 4ciC?3 + G?4 



(30) 
(31) 
(32) 



etc. The requirement that connected 3-dipole correlations vanish leads to 
the equation d 3 = c| — 1 — 3(c 2 — 1) =0 which fixes c 2 to be either 1 (no 
correlations at all) or c 2 = 2 (the solution found in [TU] ). We note that 
from the general formula (j2Sj) that solution (for (jVi) and (M1N2)) can be 
extended to a solution of the full dipole hierarchy but with some nonvanishing 
higher-dipole correlations (like ^ 7^ 0). 

5 Summary 

In this paper we studied the JIMWLK evolution equations in the approxima- 
tion of keeping only dipole-like color contractions and neglecting all higher 
multipole operators. In this approximation we have derived an equation for 
the evolution of the dipole probability distribution. No further assumptions 
like neglecting correlations were made. The hamiltonian turned out to be first 
order in functional derivatives in contrast to the full JIMWLK evolution. It 
would be interesting to study further properties of this hamiltonian. 

In the second part of the paper we derived a one-parameter family of 
exact solutions to the dipole evolution hierarchy. They are characterized by 
an enhancement (or depletion) factor in multidipole correlation functions and 
the single dipole densities are governed by a modified BK equation. 

It seems that it would be very interesting to study further the proper- 
ties of saturation when restricted to the dipole sector since this is a much 
simpler system than the full JIMWLK one. In particular one could explore 
the probabilisitic interpretation of the dipole evolution equations (fTTj). the 
structure of (attractive) fixed points and the question whether multipole op- 
erators present in the full JIMWLK framework lead to qualitative or just 
quantitative changes in the physics of saturation. 
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